Identification of RdRp inhibitors against SARS-CoV-2 through E-pharmacophore-based virtual screening, molecular docking and MD simulations approaches.

The outbreak of novel Coronavirus, an enduring pandemic declared by WHO, has consequences to an alarming ongoing public health menace which has already claimed several million human lives. In addition to numerous vaccinations and medications for mild to moderate COVID-19 infection, lack of promising medication or therapeutic pharmaceuticals remains a serious concern to counter the ongoing coronavirus infections and to hinder its dreadful spread. Global health emergencies have called for urgency for potential drug discovery and time is the biggest constraint apart from the financial and human resources required for the high throughput drug screening. However, computational screening or in-silico approaches appeared to be an effective and faster approach to discover potential molecules without sacrificing the model animals. Accumulated shreds of evidence on computational studies against viral diseases have revealed significance of in-silico drug discovery approaches especially in the time of urgency. The central role of RdRp in SARS-CoV-2 replication makes it promising drug target to curtain on going infection and its spread. The present study aimed to employ E-pharmacophore-based virtual screening to reveal potent inhibitors of RdRp as potential leads to block the viral replication. An energy-optimised pharmacophore model was generated to screen the Enamine REAL DataBase (RDB). Then, ADME/T profiles were determined to validate the pharmacokinetics and pharmacodynamics properties of the hit compounds. Moreover, High Throughput Virtual Screening (HTVS) and molecular docking (SP & XP) were employed to screen the top hits from pharmacophore-based virtual screening and ADME/T screen. The binding free energies of the top hits were calculated by conducting MM-GBSA analysis followed by MD simulations to determine the stability of molecular interactions between top hits and RdRp protein. These virtual investigations revealed six compounds having binding free energies of -57.498, -45.776, -46.248, -35.67, -25.15 and -24.90 kcal/mol respectively as calculated by the MM-GBSA method. The MD simulation studies confirmed the stability of protein ligand complexes, hence, indicating as potent RdRp inhibitors and are promising candidate drugs to be further validated and translated into clinics in future.

Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7.

New variants of SARS-CoV-2 are being reported worldwide. The World Health Organization has reported Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2) and Omicron (B.1.1.529) as the variants of concern. There are speculations that the variants might evade the host immune responses induced by currently available vaccines and develop resistance to drugs under consideration. The first step of viral infection in COVID-19 occurs through the interaction of the spike protein's receptor-binding domain (RBD) with the peptidase domain of the human ACE-2 (hACE-2) receptor. This study aims to get a molecular-level understanding of the mechanism behind the increased infection rate in the alpha variant. We have computationally studied the spike protein interaction in both the wild-type and B.1.1.7 variant with the hACE-2 receptor using molecular dynamics and MM-GBSA based binding free energy calculations. The binding free energy difference shows that the mutant variant of the spike protein has increased binding affinity for the hACE-2 receptor (i.e. ΔG(N501Y,A570D) is in the range -7.2 to -7.6 kcal mol-1) and the results were validated using Density functional theory. We demonstrate that with the use of state-of-the-art computational approaches, we can, in advance, predict the virulent nature of variants of SARS-CoV-2 and alert the world healthcare system.

Meteorological Factors and the COVID-19 Pandemic: The Backdrop of Pakistan

Coronavirus victims have been confirmed all around the world and millions of people are being put into self-isolation. In this backdrop, a superior appreciation of the effective parameters in epidemic spreading can cause a cogent assessment toward COVID-19. In this vein, the consequences of weather indicators on the spread of COVID-19 can play an instrumental role in the current coronavirus situation enveloping the world. These elements entail time, maximum and minimum temperature, humidity, wind speed, and rainfall. By such an incorporation, their consequent effects on coronavirus in Pakistan are explored. In the current study, principal elements are considered including the number of infected patients with coronavirus in Pakistan. The autoregressive distribution lag (ARDL) approach is used to analyze the effects and relationships of variables with the COVID-19 expansion rate extracting data from April 1, 2020 to April 30, 2021. The results revealed that maximum and minimum temperature, humidity, wind speed, and rainfall had a significant positive correlation with total and confirmed cases of COVID-19. Lastly, this brief communication attempts to clarify the outbreak of coronavirus in the region.

Antiviral activities of flavonoids.

Flavonoids are natural phytochemicals known for their antiviral activity. The flavonoids acts at different stages of viral infection, such as viral entrance, replication and translation of proteins. Viruses cause various diseases such as SARS, Hepatitis, AIDS, Flu, Herpes, etc. These, and many more viral diseases, are prevalent in the world, and some (i.e. SARS-CoV-2) are causing global chaos. Despite much struggle, effective treatments for these viral diseases are not available. The flavonoid class of phytochemicals has a vast number of medicinally active compounds, many of which are studied for their potential antiviral activity against different DNA and RNA viruses. Here, we reviewed many flavonoids that showed antiviral activities in different testing environments such as in vitro, in vivo (mice model) and in silico. Some flavonoids had stronger inhibitory activities, showed no toxicity & the cell proliferation at the tested doses are not affected. Some of the flavonoids used in the in vivo studies also protected the tested mice prophylactically from lethal doses of virus, and effectively prevented viral infection. The glycosides of some of the flavonoids increased the solubility of some flavonoids, and therefore showed increased antiviral activity as compared to the non-glycoside form of that flavonoid. These phytochemicals are active against different disease-causing viruses, and inhibited the viruses by targeting the viral infections at multiple stages. Some of the flavonoids showed more potent antiviral activity than the market available drugs used to treat viral infections.

Computational Investigation of Increased Virulence and Pathogenesis of SARS-CoV-2 Lineage B.1.1.7 (preprint)

New variants of SARS-CoV-2 are being reported worldwide. More specifically, the variants reported in South Africa (501Y.V2) and United Kingdom (B.1.1.7) were found to be more contagious than the wild type. There are also speculations that the variants might evade the host immune responses induced by currently available vaccines and develop resistance to drugs under consideration. The first step of viral infection in COVID-19, occurs through the interaction of receptor binding domain (RBD) of the spike protein with peptidase domain of the human ACE-2 (hACE-2) receptor. So, possibly the mutations in the RBD domain of spike protein in the new variants could modulate the protein-protein interaction with hACE-2 receptor leading to the increased virulence. In this study, we aim to get molecular level understanding into the mechanism behind the increased infection rate due to such mutations in these variants. We have computationally studied the interaction of the spike protein in both wild-type and B.1.1.7 variant with hACE-2 receptor using combined molecular dynamics and binding free energy calculations using molecular mechanics-Generalized Born surface area (MM-GBSA) approach. The binding free energies computed using configurations from minimization run and low temperature simulation show that mutant variant of spike protein has increased binding affinity for hACE-2 receptor (i.e. ΔΔG(N501Y,A570D) is in the range −20.4 to −21.4 kcal/mol)The residue-wise decomposition analysis and intermolecular hydrogen bond analysis evidenced that the N501Y mutation has increased interaction between RBD of spike protein with ACE-2 receptor. We have also carried out calculations using density functional theory and the results evidenced the increased interaction between three pairs of residues (TYR449 (spike)-ASP38 (ACE-2), TYR453-HIE34 and TYR501-LYS353) in the variant that could be attributed to its increased virulence. The free energies of wild-type and mutant variants of the spike protein computed from MM-GBSA approach suggests that latter variant is stable by about −10.4 kcal/mol when compared to wild type suggesting that it will be retained in the evolution due to increased stability. We demonstrate that with the use of the state-of-the art of computational approaches, we can in advance predict the more virulent nature of variants of SARS-CoV-2 and alert the world health-care system.